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Quantum transport through a double Aharonov-Bohm-interferometer in the presence of Andreev 
reflection is investigated in terms of the nonequilibrium Green function method with which the 
reflection current is obtained. Tunable Andreev reflection probabilities depending on the interdot 
coupling strength and magnetic flux as well are analysised in detail. It is found that the oscillation 
period of the reflection probability with respect to the magnetic flux for the double interferometer 
depends linearly on the ratio of two parts magnetic fluxes n, i.e. 2(n + 1)tt, while that of a 
single interferometer is 2-k. The coupling strength not only affects the height and the linewidth of 
Andreev reflection current peaks vs gate votage but also shifts the peak positions. It is furthermore 
demonstrated that the Andreev reflection current peaks can be tuned by the magnetic fluxes. 
PACS: 73.23.-b; 74.45. +c; 73.63.Kv 

I. INTRODUCTION 

Quantum dots (QDs) containing discrete energy levels which can be tuned easily in experiment have been attracting 
considerable attention continuously, and the study of quantum transport properties through QDs has become an 
interesting field which serves as a testing ground for the physical phenomena, such as Kondo [1,2] and Fano effects 
[3] . Recent developments on the study include Coulomb on-site repulsion [4] , spin-flip effect caused by the magnetic 
leads [5] or rotating magnetic field [6], and phonon-assisted tunneling in the molecular QD [7] with time-dependent 
field or ac gate voltage [8,9]. 

Electron transports coherently through an open QD can be investigated in the standard procedure by means 
of nonequilibrium Green function (NGF) [8], the rate equation [10] and some other methods [11]. The current 
characteristic properties change dramatically if the QD is connected to additional devices, so various kinds of hybrid 
QD systems have been studied in the literature. Two aspects mainly involved in those hybrid systems are the 
ferromagnetic and superconductor properties which evoke the spin degree of freedom of electrons. The F-QD-F system, 
i.e., QD coupled with two ferromagnetic leads, has been studied as spin-valve device [5,12]. In the superconductor 
case, the Andreev reflection, conserving the momentum of the reflected particle but not the charge, plays a dominant 
role at the interface between a normal metal and a superconductor. Using multichannel S-matrix method, Beenakker 
demonstrated the resonant Andreev tunneling in the normal metal-normal QD-superconductor (N-NQD-S) hybrid 
system in the zero-bias limit [13]; later Claughton et al. extended this theory to the finite-bias case [14]. Sun et al. 
explored the resonant Andreev reflection in this system with multienergy levels in normal QD and found that some 
extra peaks originating from the Andreev reflection superimposed on the conventional current plateaus [15]. And 
they also investigated the photon-assisted Andreev and normal tunnels in this system. With an applied time-varying 
external field, it was observed that a series of photon-assisted Andreev tunneling peaks appeared in the curve of 
the average current (I) vs the gate voltage with the negative peaks on the left-hand side and positive peaks on the 
right-hand side of the original resonant peak in the absence of the external field [16]. Another interesting hybrid 
system consists of ferromagnet and superconductor leads coupled with a normal QD (F-QD-S). In Ref. [17], Zhu et 
, al. provided an efficient method for introducing a spin degree of freedom in this system, and also derived the general 
current formula of the Andreev reflection in a quantum dot connected with two ferromagnetic and one superconductor 
• i-H ! leads (Fi, F 2 -QD-S) [18]. Considering both the Coulomb interaction and the multienergy levels in QD, Feng and 
Xiong studied the transport properties in the F-QD-S system [19]. Very recently, Cao et al. considered the spin-flip 
effect in QD of a F-QD-S system and found that the single-resonance peak of the Andreev reflection conductance 
versus the gate voltage would spill into a double-resonance peak with the increasing of the spin-flip scattering strength 
[20]. Thus far, most of studies relating to Andreev reflection have been focused on a single-QD system. 

On the other hand, the electron transport through a double-quantum-dot (DQD), especially the parallel-coupled 
double-quantum-dot (PDQD), is an interesting subject which may involve the quantum phase interference. An open 
PDQD threaded by a magnetic flux is regarded as an Aharonov-Bohm (AB) interferometer in which the coherence of 
the electron remains partly even in the presence of Coulomb repulsion [1,21]. Tunable coupling strength between two 
parallel dots makes the transport properties more complex and interesting [22,23]. The tunable Andreev reflection in 
an AB interferometer has been investigated recently by Peng et al. [25]. 
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In this paper we shall study the electron transport through a N-PDQD-S structure, i.e., a parallel-coupled double- 
quantum dot that is connected with two leads, one is a normal metal and the other is a superconductor. The device 
threaded by two magnetic fluxes (see Fig.l) with the interdot coupling can be regarded as a double AB interferometer. 
By using the NGF method, the current Jl and the probability of the Andreev reflection are derived. 

The paper is organized as follows: In Sec. II, we present the model Hamiltonian and derive the general formula of 
the current and the probability of the Andreev reflection with NGF method. In Sec. Ill, the properties of the current 
due to Andreev reflection depending on bias voltage, magnetic flux, gate voltage, and interdot coupling strength are 
studied in detail, and several dufferent effects are predicted. Finally, a brief summary is given in Sec. IV. 



II. MODEL HAMILTONIAN AND FORMULATION 



The Hamiltonian of the system schematically shown in Fig. 1 can be decomposed into four parts of the following 
form: 

H = H L + H R + H D + H T , (1) 

where 

H L =^2(£k- eV)a\ a a k a, 

kcr 

H R = Y £ P b U b P° + ^( Ab Pi b -P^ + h - c -)> 

pa p 

H D = e t dld ia +Y / (^^4 a d la + h.c), (2) 

cr, i=l,2 a 

H T = ]T (Lie^aljic + h.c.) 

fc<T,l=l,2 

+ W* R bl a d la + h.c.). 

pcr^i— 1,2 

Hl describes the left normal metal lead, atcr is the corresponding annihilation (creation) operator, and V is 

the bias voltage. We assume the chemical potential of right lead — [16]. Hr describes the right superconductor 
lead with the energy gap A. Hd models the coupled double quantum dots, is the interdot coupling strength, and 
<f>i2 is the corresponding phase shift induced by the magnetic flux. Here we consider the single energy level in each 
quantum dot. H T is the tunneling part of Hamiltonian, Li (Ri), independent of momentum k (p), is the hopping 
strength between ith quantum dot [i. e. the quantum dot with energy level Si (i = 1,2)] and left (right) lead, and 
$li is the phase shift of the electron tunneling from the left lead to first quantum dot, where 

$ 

$L1 = &1R = —r 

4 

where $ = 27r(ai + a 2 ) with $^ = ai$o (« = 1,2) being the magnetic flux in the quantum unit $ = hc/e and 
$12 = -$21 = 2n(ai - a 2 ) . 

The current from the left lead to a quantum dot can be written as [4,8] 



J ij = l,2 kcr 

xS(e-e k ) [°° dt^G^t^fde + eV) 

J —OO 



+e(t'-t 1 )G< <J (t,t 1 )}e-^ t '\ (3) 

where we define G r ij(7 (t,t') = -i6(t - t') ({d icr (t), and G< JrJ {t,t') = i (d) a {t')d ici {t)^. Because of the right 

superconductor lead, it is convenient to introduce 4x4 matrix representation in which G r (t, t 1 ) and G < (t, t 1 ) have the 
forms, 

G r {tX)^-i8{t-t')({^(t)^\t')}), (4) 
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G<(*,* / ) = i<* t (* / )*(*)>, 
where ^ — (d\^, d\i, d\^, d 2 j.). So Eq. (3) can be rewritten as 

JL T (t) = -f Im /fj/* A'e^ t - t '){r L [G r (t,i')/L(e + ^) + G<(t,i')]}ii+33, 



(5) 



(6) 



where 



rf 
o 



/ LiLJ I^e 4 * \ 

L X L\ LiL^e^t 

LiL^e"^ L 2 U 2 

J L 2 i 1 e J * L 2 L* 2 J 

v/ffffe 4 * \ 



V 



L e l ~ 







rf 
o 



1 2 



e 2 





1 2 



(7) 



/ 



with rf = 2nJ2k — £k)L\L\, T 2 = 2nJ2k — £k)L 2 L 2 , here we consider only a spin-up current Jj,|(t) [ </lj,(£) 
is easily obtained from Jlt(^) by exchanging the spin index]. 
To calculate G r (i,i'), we start from the Dyson equation, 



~(t,tf) = g r (t,t') + J di 1 dt 2 G r (i,ii)S r (i 1 ,i 2 )g r (i 2 ,i'), 



(8) 



in which g r (t, t') is the retarded Green function for the isolated quantum dots, and can be easily obtained as, 



/ g — »£l(t — *') 



? r (i,t') = -i6(t-t') 



V 





fei(t-t') 









e -ie 2 (t-t') q 

o e *e2(t-t') y 



(9) 



Self-energy £ r (£i,i 2 ) consists of three parts, such that, 

E r (ti, t 2 ) = S 12 + S£(ti,t 2 ) + Ek(ti,t 2 ), 



(10) 



where E£(ii,f 2 ) results from the tunneling coupling between QD and the normal metal lead, l¥ R {ti,t 2 ) between QD 
and superconductor lead, and Si 2 is from interdot coupling contribution. Under the wide-bandwidth approximation 
[8,16,20], the self-energy becomes 
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and E£(ii,t 2 ) is written as 
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The self-energy from the coupling between QD and right superconductor lead is given by [16], 
SJj(ti,t 2 ) 
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and the retarded Green function of the isolated superconductor lead glp{t,t') is read as 
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with p 1 ^ being the density states of the right lead in a normal metal state. We define rf = 2ir \ Ri | 2 p 1 ^. By substituting 
Eqs. (9)-(13) into Eq. (8), the Green function G r (t, t') can be derived. It is, however, convenient to transform Eq. 
(8) into the energy representation by making use of the Fourier transformation, 



G r (E,E')= J dtdt'G r (t,t')e lEt e- lE,t ' . 

Then solving Eq. (8) in the energy representation, we get (see Appendix A) 

G^ j (E,E') = 2nG r ij (E')5(E-E'), 



(17) 
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here, i,j = 1,2,3,4. G\j{E,E') arc the elements of the matrix G r (E,E'), and the detailed calculation of G\j{E') is 
given in Appendix A. The Green function G r (t,t') is obtained after making an inverse Fourier transformation as 
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In the following, we derive the smaller Green function using the Keldysh equation 

G < {t,t') = J dt 1 dt 2 G r (t,t 1 )-E < (t 1 ,t 2 )G a {t2,t'), 
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Pr(s) is the ratio of superconducting density of states to the normal density of states [15]. According to Eq. (20), 
G < (t,t) can be obtained by substituting Eqs. (21)-(23) and the Green functions G r (ti,t 2 ) , G a (ti,t 2 ) into Eq. (20) 
with the results of matrix elements Gf x (t, t),Gf 3 (t, t), G 31 (t, t), and G 33 (t, t) shown in Appendix B. 

The total current, i.e., the sum of both spin-up and spin-down parts, is found by substituting the retarded Green 
function G\At,t') and the smaller Green function GfAt,t) shown in Appendix B into Eq. (6), and the result is 



Jl = J A + Jt 



(1) 
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Ja is seen to be the current due to the Andreev reflection contribution and the probability of the Andreev reflection 
T AR is given by 



+rfrf ( G[ 4 (s) + G r 32 (e) ) 
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As explained in Refs. [15,24], the term (rf) 2 G\ 2 {e) [(r^) 2 G 34 (e) ] describes an electron with spin-up tunnelling 
from the left lead through quantum dot 1 (2) to the superconductor lead and reflecting a hole back to the left lead 
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through the quantum dot 1 (2). While the term TfTf' G r u (e) \T{Y% G r 32 (e) ] is for the spin-up electron tunneling 

from left lead through dot 1 (2) to the superconductor lead and a hole reflecting back to the left lead through the 
dot 2 (1). There exist other channels in which the tunnel electron and reflecting hole go through both two dots. 
The remaining part in Eq. (26) describes the interference effect of all possible multichannels of the device. The 
transmission probability from left to right leads reads 



r M = rf[rf G^e) + rf G[ 2 (s) 
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+2Re[-T^e- i i^-G r u (s)G{* 2 (6) 



containing many terms which are understood. We are mainly interested in the physics of the Andreev reflections. As 
a matter of fact only the Andreev reflection process contributes to the current (i. e., Jt — 0) at zero temperature, if 
\eV\ < A. In the following section, we analyze time evolution of the probability of the Andreev reflection Tar and the 
Andreev reflection current J a as well depending on dc bias voltage, gate voltage, magnetic flux, and interdot coupling 
strength. 



III. NUMERICAL RESULTS 



First of all we transform the coupled-QD system into an effective decoupled one with bonding and antibonding 
dressed states of the QD molecule using the following transformation [23], 

( f+ \ _ ( - cos^e 4 * 12 - sin/3 
\ f- ) ~\ - sin p cos /3e- 1 * 12 

where, /_ and /+ are the annihilation operators of the bonding and antibonding dressed states for the QD molecule. 
And (3 is defined as l/2tan _1 [2ft/(£i-£ 2 )]. Thus the Hamiltonian for coupled double quantum dot Hp is diagonalizcd 
as 

H D = E + flf + + E_flf_, (29) 

where 

E± = i[£i + £ 2 ± ^(ei - £ 2 ) 2 + 4ft 2 ]. (30) 

Considering the tunneling Hamiltonian between double QD and the left lead, we obtain the linewidth and 
corresponding to the antibonding and bonding states coupled to the left lead, respectively, as 

rf = rf cos 2 p + sin 2 p + yrfr^ sin(2/3) cos (| + $ 12 ), (3i) 
= sin 2 p + cos 2 p - yrfrjsm(2/?) cos (* + $ 12 ). (32) 

In the numerical analysis we set A = 1 as the energy unit. The energy dependence of reflection probability Tar 
on the interdot coupling strength is shown in Fig. 2 for the symmetry [with = Tf = T 2 = T R = 0.02 , Fig. 2(a)] 
and asymmetry [with rf = T% = 0.08, rf = T R = 0.02 ,Fig.2(b) and = rf = 0.08, T 2 = T R = 0.02 ,Fig.2(c)j 
cases, respectively. Other parameters are chosen as e\ = e 2 = 0, and <f> = <f>i 2 = 0. It is seen in Fig. 2(a) that if 
the two QDs arc completely decoupled (i.e., O = 0) only a single peak appears (solid line), and when ft ^ , two 
peaks emerges at e = ±Q due to the existence of two energy levels. The maximum of peaks is getting lower when the 
coupling constant ft is increasing. For the asymmetry case Fig. 2(6) shows that there are no significant two peaks. 
Fig. 2(c) exhibits the complex spectra due to the possible overlap of the broad levels similar to the interpretation of 
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Fano effect in a parallel coupled double QD connected to two normal metal leads system [23]. And when fl is large 
enough such that the broad two energy levels are completely separated, the two single peaks remain at e — E + and 

It is interesting to see the effect of the magnetic flux on the Andreev reflection and Fig. 3 presents the dependence 
of the Andreev reflection probability on the magnetic flux $ with ct\ — a 2 and therefore $12 = 0. The height of 
peaks at e = ±f2 can reach 2 when the magnetic flux is nonzero, and magnetic flux also has an influence on the 
width of the peaks, i. e., on the lincwidth function T±, here we choose $ = tt/3 for a solid line, and 27r/3 for a dotted 
line, the height of peaks are not changed. Electron phase coherence is shown in Fig. 4. When ft — 0, i. e., there is 
no the interdot coupling, the oscillation period of the Andreev reflection probability vs magnetic flux is 27r$o (solid 
line), while U / 0, the period changes to 47r<f> when a\ = a 2 [see Fig. 4(a)]. Here, we choose the different energy 
eigenvalues such that e = 0.11, since at those values the probability of the Andreev reflection reaches maximum. To 
sec the reason of the period changing, we give a general result of the oscillation period when ai/a 2 = n. The channel 
paths contain two parts due to the Andreev reflection in the hybrid system. One is the incident electron from the left 
lead to the superconductor in Fig . 5a, the other is the reflecting hole from the superconductor to the left lead in Fig. 
56. The phase shift of electron in the case (3) in Fig. 5a is ($1 — $2)/2, and the corresponding phase shift of hole in 
the case (1) in Fig. 56 is ($1 + <& 2 )/2, so the total shift of the channel (3) for the elcetron and (1) for the hole is $1. 
The interference effect corresponds to all the channel paths reads [25] 

A + Aie 4 * 1 + Aae-'* 1 + + A 4 e-^ + A 5 e t{ ^+^ 2 

= B + Bx cos $1 + B 2 cos $ 2 + B 3 cos($i + $ 2 ) 
+B 4 cos($i - $ 2 ) + B 5 cos 2$i + B 6 cos 2<J> 2 
+B 7 cos 2($i + $ 2 ) + B$ cos(2$i + $ 2 ) + B 9 cos($i + 2$ 2 ) 
+B W cos 2($i - $ 2 ) + Bn cos(2$i - $ 2 ) + B 12 cos($i - 2$ 2 ) 

(33) 

where parameters A are the amplitudes of channel paths. With $1 = $n/(n + 1),<&2 = $/(n + 1), Eq. (33) shows 
the oscillation periods might have ir, 2tt, 2ir(n + l)/n, 2ir(n + 1), 2n(n + l)/2n, 2tt(ti + l)/2, 2n(n + l)/(n — 1), 
2?r(n + l)/(n - 2), 2?r(n + l)/(n + 2), 27r(n + l)/(2n + 1), and 2?r(n + l)/(2n - 1), suggesting that the oscillation 
period of reflect probability should be 2{n + The linear depedence of the ratio n is further confirmed in Fig. 
4(b), it shows that the oscillation period is 6ir for n = 2 (solid line), 87r for n = 3 (dotted line), and 107r for n = 4 
(dashed line). The linear manner of the osicllation period is the same as the one in Ref. 22, however, the mechanism 
is different: there the transmitted current is analyzed. 

Fig. 6 shows current-gate-voltage curve (Ja vs V g ) for the case that E\ — 0.1, e 2 — 0.3, $ = $12 = 0, and the 
bias voltage V — 0.4. We assume that the energy levels in both dots are connected to the same voltage through the 
capacitive coupling. When O = (solid line), there exist two smaller peaks at V g = 0.1, and 0.3 located symmetrically 
beside the larger peak at V g — 0.2 with the height that is two times higher than the height of the smaller peaks. The 
interpretation of the phenomena is similar to that in Ref. [15] such that when V g = 0.1 the energy level £\ matches 
with chemical potential ur = in the right superconductor lead and the Andreev reflection takes place at this level. 
When V g — 0.2 two energy levels have a symmetric localization about ur ( e 2 is 0.1 above ur while z\ is 0.1 below 
ur ) and both £\ and e 2 are below the bias voltage, so the Andreev reflection that involves those two levels leads to 
the height of the peaks that are two times higher than that containing only one energy level. The interpretation is 
the same for the gate voltage V g = 0.3 where the energy level e 2 matches with ur. Some new properties are observed 
for the coupled case f2 7^ 0: (1) the location of the middle peak does not change, while the positions of the left and 
right peaks shift towards the left and right sides, respectively, with the increasing coupling strength Q. (2) the heights 
of the left and the middle peaks are suppressed significantly; meantime, the height of the right peak is enhanced. 
(3) the widths of the left and middle peaks decrease while the width of the right peak widened. These interesting 
phenomena can be explained with Eq. (30)- (32). The energy level of dressed state E + gets larger and becomes 
smaller with increasing coupling strength £1 , and at the end E+ is higher than energy level of bared state e 2 , when 
is smaller than s\. So when the gate voltage takes action on the dressed states, the left peak runs towards the 
left while the right peak goes opposite, however, the locations of the two peaks remain symmetric about the middle 
peak. We assume £L = e\ — a — 0.1 — a, E + = e 2 + a = 0.3 + a, where a is introduced due to the coupling strength 
fl. So when the gate voltage becomes V g = 0.1 — a, the left peak emerges and with increasing V g (up to 0.2), both two 
dressed states levels locate symmetrically about ur. Thus, the location of middle peak does not change at V g = 0.2, 
no matter how the parameter a(fl) varies. The linewidth of dressed states, , increases with the increasing of f2, 
while decreases under the parameters considered here, and hence, the right peak gets wider and left peak becomes 
narrower with the increasing of fi. The linewidths T± are directly connected with the coupling strengths between 
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QD and the leads. The widening of the right-side peak means the increasing of coupling strength between dressed 
state E + and the leads. Therefore, the height of the right-side peak is enhanced, and opposite to this, the height 
of the left-side peak is suppressed. The suppression of the middle-peak height suggests that the Andreev resonance 
involving two energy levels gets maximum when two energy levels are symmetric. 

Figures. 7 — 9 show a tunable Andreev reflection current depending on the magnetic flux when fl ^ 0. In Fig. 
7, we fix $ = 7r/2, and tune the value of $12 for different curves. One can see that the height of the right peak is 
enhanced and the height of the left one is suppressed while the height of the middle peak does not change with the 
increase of $12. It is interesting to notice that the left peaks has essentially no changes when $12 = (solid line) 
and 7r/2 (dotted line), while the right peak changes conspicuously with respect to the case without the magnetic flux. 
Figure 8 has the opposite effect compared with Fig. 7, namely, the height of the left peak is suppressed while the 
height of the right one is enhanced and the middle peak does not vary with $ 12 . We fix $ 12 = ir/5, and vary values 
of magnetic flux $ in Fig. 9, it is shown that the height of three peaks can be tuned simultaneously. 



We have studied the tunable Andreev reflection in a double AB interferometer in terms of the noncquilibrium Green 
function method and observed several features of the reflection current. It is found that the oscillation period of the 
Andreev reflection probability with the magnetic flux is 2ir when interdot coupling vanishes ( = 0, in this case the 
system reduces to a single interferometer), while it is 2(n + l)ir for our double interferometer ( SI 7^ 0) where n is the 
ratio of two parts magnetic fluxes, i.e. n = $i/<f>2. The Andreev reflection current peaks cannot only be tuned by 
coupling strength fl, but also by magnetic flux. 



In this appendix, we present the derivation of the Green function of Eq. (18) in detail. From Eqs. (8)-(13) and 
(17) we have the following formulas for the matrix elements of Green function G r (t, t') in the energy representation: 



IV. CONCLUSIONS 
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VI. APPENDIX A 
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4 1 = l(rf + rf , ^ 

2 V 1 1 ^{ E 'Y - A 2 



1^1 

V(^') 2 -A 2 



2 V(^') 2 -A 2 ' 



and 



si = fie'* 12 - ^rf r£ e -** - * 



|E'| 



2 ^(S') 2 - A 2 



rfr 2 Vf, 



x 2 = -fi*e- J * 12 - ^rfT^e^ - 



2 V(^') 2 - A 2 

2V 1 2 2 ^(S') 2 - A 2 V 

l#l 



X4 



2V" 1 " 2 " 2 V /(^')2_A 2 
Solving those coupled equations, the results are obtained as 

GI 1 (£,^) = 27rG; 1 (^)5(£;-^) ) 
GfI 2 (£;,^) = 27rGI 2 (^)5(f;-^) ) 
G[ 3 (£;,^) = 27rG[3(^)5(£;-^) ) 
G r 14 (E, E 1 ) = 2nG r 14 (E')S(E - E'), 



rfife**. 



where 



GTi(^) 



M(E') 



g^ 1 (E')+A 1 -Bv/ffTt 



M(E') 



-BTfe'^ -By/Tfff -x 2 

-x 3 g^\E')+A 2 -BT§e-^ 
-By/fM ST^i g^\E')+A 2 



G r 13 (E>) 



M{E>) 



-X3 -s^rfrf 



-X2 

-BT^e~ 



hcrc, 



-£ 4 



-BTffe 1 - 



-X3 



M(E') = 
With the same procedure we have 



g r ^\E') + A 1 -flrfe* 



-x 4 -srfe'f ^^^o + A 2 



g r n\E')G r 31 {E, E>) = -A x G r 31 (E, E') + BT?e^G r 32 (E, E') 
Xl G r 33 (E,E>) + B^T^G r 34 (E,E'). 



9 r 2 2\E')G r 32 (E,E') = Br?e-*-G r 31 (E,E>) - A 1 G r 32 {E,E') 

+B^T^G r 33 {E,E')+x 2 G r 34 {E 1 E l ) 1 

9 r 33 \E')G r 33 (E,E') = 2irS(E-E')+x 3 G r 31 (E,E') 

+B^T^G r 32 (E, E') - A 2 G r 33 (E, E') 
+BT^e- i iG r 3i (E,E'), 



So the final results are 



where 



gl^(E')G r 34 (E,E') = B\JvfYfG 31 (E, E') + x 4 G 32 (E, E') 
+BT^iG r 33 (E,E') - A 2 G r 44 (E,E'). 



G r 31 (E,E') = 2nG 31 (E')5(E-E') 
G r 32 (E,E') = 2irG r 32 (E')8(E - E') 
G r 33 (E,E')=2nGl 3 (E')S(E-E') 
G r 34 (E,E') = 2nG r 34 (E')5(E-E') 



G r 31 (E>) = 



M(E') 



-srf 



e 2 



-Xl 



-b^tM 

g r 22 1 {E') + A l -Bv/Tfrf 

-x 4 -BT^i g r 44 1 (E') + A 2 



M(E') 



g^iE^+A, -x x -S^ffTf 
-BY fe-^ -B^/fftf -x 2 
-B^ffff -BT§e^ g r 44 1 {E') + A 2 



G r 33 (E') 



M(E') 



g r ^ 1 {E')+A 1 -Blfe*? 



OTfe-^ g r 22 1 (E') + A 1 



-x 2 

7 1 / rpl\ 



-x 4 g r 44 i (E')+A 2 



G 34 {E') 



M(E') 



g r v [ 1 {E')+A 1 -BT*e l i -xi 
-Srfe"** g r 22 1 {E') + A l -B^fftJ 
-B v / Tff 2 1 -x 4 -BT§e l i 



VII. APPENDIX B 

The complete expression of matrix elements of the Green function G < (t, t) is 
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Gf 1 (t 1 t)= l J ^.{f L{ s + eV)r^\G r n (e)\ 2 + f L (e-eV)Ti 

+f L (e + eV)r% \G r 13 (e)\ 2 + f L (e - eV)T% \c r 14 (e) 

+2^TfT^ Re[f L (e + eV)e l i G r n (e)Gr 3 (e) 
+f L (e-eV)e-^Gl 2 (e)G\l(e)]} 

+l / £/«(*( £ ){ r ?| 5 n( £ )r+ r ? l^iaOO 



i 



+if 



Gi 3 (e) 



A 



+2Re[-rfe-^ n G^( £ )G^( £ ) 



+v /rfrf e -^G^( £ )GS( £ ) + v^fr? e^GJ 2 ( £ )£?S( e ) 



-If e't ^GI 3 ( e )Gi:( £ )) - VTfrf-G^( e )Gn(e) 

-V / rfr? A^ 2 ( £ )G^( £ )]}, 



Gf 3 {t,t) = i f ^[J L (e + eV)Y L 1 Gl 1 (e)GUe) + h(e + eV)^^e-^G[ 3 (e)GUe) 



+.f L (e - eV)TfG[ 2 (e)G a 23 (e) + f L (e - eV)^T^e^G r 14 (e)G a 23 (e) 
+f L (e + eV)^T^e l iG r n (s)G a 33 (e) + f L (e + eV)T^G r 13 (e)G a 33 (e) 



+f L (e - eV)JT^e-^Gl 2 (e)Gt 3 (e) + f L (e - eV)Y^G r u (e)Gl 3 {e)] 



A 



- v / rfrf A^ 4 ( £ )]g? 3 ( £ ) + [-rfe-^^G^( £ ) + rf gj 2 ( £ ) 



;A 



-V r f lf^ 3 ( £ ) + V r ? r«e-^G^ 4 ( £ )]G5 3 ( £ ) 

+[V / rfrf e-'T g^( £ ) - ^/rfrf ^G[ 2 (s) + i?gj 3 ( £ 

A 



T*e-^-GUe)]G a 33 {e) + [-Jt^-G^s) 



+ 



v / rfrf e^G' 12 ( £ ) - rf e'f ^g^ 3 ( £ ) + r^ 4 ( £ )]G^ 3 ( £ )}, 



Gf 1 (t,t) = i / J[/L( £ + e y)rfGS 1 ( £ )G? 1 ( £ ) + / i ( £ + e y)v / r7^e-^G5 3 ( £ )G? 1 ( £ ) 



+/ L ( £ - eF)rf Gg 2 ( £ )G^( £ ) + f L (e - eV)^T L 1 T^-G r 3A {e)G a 21 {e) 
+f L (s + eV)Jr^e i iG r 31 (e)G a 31 (e) + f L (e + eV)Y L 2 G r 33 {e)G a 31 {e) 



+f L (e - eV)y/r^e-^G r 32 (e)G a il {e) + f L (e - eV)T^G r 34 (e)G a 41 (s)} 

-V^wA^^^G^^ + l-rfe-^^G^^ + rfG^^) 
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+ [jT*r*e- i *G r 31 (e) - Jr? if^G^e) + if G r 33 (e) 

+Vrfrf e^G5 2 ( £ ) - e ^pG5 3 ( e ) + rf g^g^)}, 



G&(t,t) 



'i(e + eV)rf G^( £ ) +f L ( E -eV)T[ G r 32 (e) 



+ f L (e-eV)T% G r 34 (e) 



+2y/T^ Re[f L (e + e^G^G^fe) 
+f L (e-eV)e-^G r 32 (e)G r 3 l(e)}} 

/dp ~ 2 ~ 

-/«( £ )^( £ ){rf g^( £ ) +rf g^ 2 ( £ ) 



+r 



G r 33 (e) 



+ rf 



+2Rc[-rfe^^G^ 2 ( £ )G^( £ ) + y/T* If e^G^( £ )G^( £ ) 
+ V / Tfrf e^G5 2 ( £ )G^I( £ ) - r«e^^G53( e )GS( e ) 



(B4) 



The results here show that the above matrix elements of G < (t, t) do not depend on time t, because for the N-QD-S 
hybrid system, the current should be independent of time for the dc bias [15]. 
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FIGURE CAPTION 

FIG. 1: Schematic diagram for the double AB interfermeter connected with normal metal and the superconductor 
leads, respectively. 

FIG. 2: The probability of the Andreev reflection Tar vs e in the unit of the energy gap A for (a) the symmetry 
case( rf = rf = rf = rf = 0.02) and asymmetry case (b) Y[ = rf = 0.08,rf = rf = 0.02, (c) rf = Tf = 0.08, 
Tf = = 0.02. The solid, dotted, dot-dashed, and dashed lines correspond to ft=0, 0.02, 0.04, and 0.08, respectively. 

FIG. 3: The magnetic flux $ dependence of the Andreev reflection probability for rf = Tf = rf = rf = 0.02, 
£1 = £2 = 0, ft = 0.2, and $12 = 0, solid line: $ = tt/3, dotted line: $ = 2tt/3. 

FIG. 4: The periodic oscillation of the Andreev reflection probability Tar with the magnetic flux <f> ( Tf = rf = 
rf = Tf = 0.02, £1 = -0.1, e 2 = 0.1), (a) oi/a 2 = 1. Solid line: ft = 0, e = 0.1, dotted line: ft = 0.05, e = 0.11, 
and dashed line: ft = 0.1, £ = 0.14. (b) ft = 0.05, £ = 0.11. Solid line: 0:1/0:2 = 2, dotted line: oi/o 2 = 3, and 4 for 
dashed line. 

FIG. 5: Paths and corresponding phase shift, (a) for the incident electron from left lead to superconductor (b) for 
the reflecting hole from superconductor to left lead. 

FIG. 6: Andreev reflection current J A vs gate voltage V g ( £1 = 0.1, £ 2 = 0.3, rf = rf = rf = Tf = 0.02, and 
$ = $12 = 0) for bias voltage V = 0.4, ft = (solid line), ft = 0.05 (dotted line), and ft = 0.1 (dashed line). 

FIG. 7: Andreev reflection current vs gate voltage with the fixed total magnetic flux that <& = tt/2 for ft = 0.05. 
$12 = (solid line), 7r/2 (dotted line), 37r/4 (dot-dashed line), n (dashed line). 

FIG. 8: Andreev reflection current vs gate voltage for $ = 37r/2. $i 2 = (solid line), 7r/4 (short-dashed line), n/2 
(dotted line), ir (long-dashed line). 

FIG. 9: Andreev reflection current vs gate voltage with $12 = 7r/5 for $ = (solid line), 27r/3 (dotted line), 57r/3 
(dashed line). 
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